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ABSTRACT 

He  use  explicit  self  similar  solutions  of  suitable  model  problems  and 
comparison  arguments  to  study  the  regularity  of  the  free  boundary  near 
t  -  0  for  the  one-dimensional  degenerate  Cauchy  problem 

vt  *  *<v>xx'  X  €  R,  t  >  0  , 

(P) 

v(x,0)  =  g(x)  , 

where  $ (v)  *  max(v,0).  The  datum  g  is  assumed  smooth  on  R\{0}  with  at 
most  exponential  growth  at  infinity  and  satisfies 

xg(x)  >0,  x  f  0,  g(0 )  -  0  , 
the  inequality  being  strict  near  x  *  0. 
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SIGNIFICANCE  AND  EXPLANATION 

Degenerate  parabolic  equations  arise  in  the  description  of  melting 
processes,  gas  dynamics  and  certain  biological  models.  The  interfaces 

ft 

corresponding  to  degeneracies  in  the  constitutive  function  usually  separate 
different  media  in  the  physical  problem. 

Problem  (P)  stated  in  the  abstract  is  related  to  nonlinear  diffusion 


equations  with  nonmonotone  constitutive  functions  as  has  been  discussed  in 


i 


J 


'  [HN1-3] y  In  this  report  we  obtain ^xplicit  self-similar  solutions  for  (P) 
corresponding  to  a  class  of  model  initial  data  and  determine  the  free  boundary 
explicitly.  The  qualitative  behavior  of  these  solutions,  in  particular  of 


their  interfaces,  is  typical  of  the  situation  in  more  general  problems.  For 

// c.  a.  u  t  !•»•>»- 

general  initial  data,  we  then  use^these  self-similar  solutions  as  comparison 


functions  to  study  the  regularity  and  the  behavior  of  the  free  boundary  for 
small  time. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  author  of  this  report. 


A  CLASS  OF  ONE-DIMENSIONAL  DEGENERATE  PARABOLIC  EQUATIONS 

John  A.  Nohel 

To  Jarosiav  Kurzweil  on  his  sixtieth  birthday 


1 .  INTRODUCTION .  In  this  paper  we  use  self  similar  solutions  to  study 
the  one-dimensional  degenerate  Cauchy  problem 

vt  =  <Mv)xx,  x  e  R,  t  >  0  , 

(P) 


v(x,0 )  *  g(x) 

for  the  constitutive  function  $ ( v)  =  max(v,0).  We  assume  that  the  initial 
data  g  are  smooth  on  R\{0}  with  at  most  exponential  growth  at  infinity  and 
satisfy 


xg(x)  >0,  x  f  0,  g(0)  ■  0  ,  (g) 

and  the  inequality  is  strict  near  x  *  0.  The  principal  issue  is  the 
regularity  of  the  free  boundary  near  t  ■  0. 

Problems  of  this  type  arise  as  convexifications  of  diffusion  equations 
with  nonmonotone  constitutive  functions  as  has  been  discussed  in  [HN1,H] •  The 
behavior  of  solutions  for  (P)  is  similar  to  those  of  the  one  phase  Stefan 
problem  in  which  g(x)  =  -1  for  x  <  0,  and  g  has  a  jump  discontinuity 
at  x  *  0.  However,  (P)  and  the  Stefan  problem  exhibit  different  regularity 
properties  of  the  free  boundary  near  t  ■  0. 

Existence  and  uniqueness  of  weak  solutions  of  (P)  can  be  proved  using 
nonlinear  semigroup  theory  [BCPa , E] .  indeed,  Benilan,  Crandall,  and  Pierre 
[BCPi]  have  obtained  optimal  existence  and  uniqueness  results  for  the  porous 
medium  equation  (i.e.  problem  (P)  with  $(v)  *  max(v|n,0),  m  >  1)  in  if1. 


However,  their  results  do  not  apply  directly  to  the  case  m  ■  1.  Details  of 


the  proof 8  in  the  present  context  are  given  in  Vizquez  [VI]  and  the  results 
(Theorems  3.1/  3.2)  are  stated  in  Section  3.  Using  standard  approximation 
arguments  and  the  comparison  method  discussed  in  [VI]  (see  Section  3)  one 
establishes  the  existence  of  a  continuous  monotone  decreasing  free  boundary 
t  s(t),  s(0)  *  0  and  v(s(t)+#t)  ■  0.  Moreover,  the  pair  (v,s)  satisfies 
the  free  boundary  problem 

vt  "  vxx'  x  >  ■<*>»  t  >  0  , 
v(x,0)  -  g(x),  x  >  0  , 

v(s(t)+,t)  ■  0  ,  (P) 

g(s(t))s'(t)  -  vx(s(t)+,t),  t  >  0  , 
s  ( 0 ;  -  0  . 

Conversely,  the  solution  v  of  (P)  extended  by  v(x,t)  ■  g(x)  for  x  <  s(t) 
is  a  weak  solution  of  (P). 

In  the  special  case  of  (P),  (g)  with  the  function  g  satisfying  the 
additional  assumption:  g'(0+),  g'(0“)  f  0,  we  have  shown  by  solving  a 
singular  integral  equation  [HN2,3]  that  the  problem  (P)  has  a  unique 


solution  ( v, s ) .  Moreover , 


$(t)  =  -x/t  +  0(t°+1/2)  (t  ♦  0+)  , 


(1.1) 


for  any  0  <  a  <  1/2,  where  <  is  the  uniquely  defined  monotone  function 
of  p  :=  g' (0+)/g' (0“)  determined  implicitly  by  the  equation 


2  3  2..  • 

-K 


e-*2/4dy 


(1.2) 


For  t  >  e  >  0  the  free  boundary  is  smooth.  It  should  be  observed  that  a 
formal  expansion  of  the  ODE  g(s(t))s'(t)  =  vx(s(t)+»t)  yields 


(p_)s(t)s' (t)  -  (p+)  +  •••  , 


so  that 


-  *  k  *  to  *  •» 


where  •••  denote  terms  of  higher  degree  powers  of  t.  However,  if  p  ■  1, 
then  (1*2)  implies  k  *  .9034  ...  f  /2.  The  reason  for  this  apparent 
inconsistency  is  that  all  derivatives  of  v  become  singular  at  (x,t)  at 
(0,0)  and  in  particular,  vx  is  discontinuous  at  (0,0).  By  contrast,  for 
the  one  phase  Stefan  problem  [FP 1 , 2 >  KM;  S]  the  solution  is  smooth  on  the 
set  {(x,t)  :  x  >  s(t),  t  e  [0,T]}. 

The  purpose  of  this  note  is  to  discuss  the  regularity  and  qualitative 
behavior  for  small  t  of  the  free  boundary  of  (P)  for  more  general  initial 


data  g.  Of  particular  interest  is  the  model  datum 


P+x\ 


x  >  0  , 


?_  |  X I Y ,  x  <  0  , 


where  p^,  y  >  0  are  given  constants,  for  which  the  integral  equation  method 
of  [HN2,3]  breaks  down  if  Y  t  1#  if  the  datum  g  is  not  piecewise 

linear.  It  will  be  shown  in  Section  2  that  (P)  with  the  model  initial  datum 
(D)  can  be  solved  explicitly  for  all  y  >  0  using  self  similar  solutions,  and 
the  free  boundary  is  determined  explicitly  for  all  t  >  0  by  proving  the 
following  result. 

Proposition  1.1.  For  the  model  datum  (D)  problem  (P)  has  the  unique  self 


similar  solution 


v( x, t )  -  tY/2*(— ),  x  >  s(t ) ,  t  >  0 

✓t 


where  ♦(•)  is  the  unique  solution  of  the  ordinary  differential  equation 


2*"<«  ♦  *f(5) 


subject  to  the  initial  conditions 


I 


<1.5) 


and  k  is  related  to  p^  via  the  condition 

lim  »  P+  . 


The  free  boundary  is  s(t)  =  -</t,  t  >  0.  Moreover/  given  p  *  p+/p.  >  0, 
k  >  0  is  the  unique  solution  of  the  equation 


(1.6) 


where  D_^_.|(»)  is  the  parabolic  cylinder  function  of  index  <-y  “  D 
[B,  p.  119]. 

It  is  readily  verified  that  (1.6)  reduces  to  (1.2)  if  y  *  1.  The 
regularity  of  the  free  boundary  for  problem  (P)  with  initial  data  (g)  more 
general  data  than  (D)  can  then  be  discussed  using  the  above  explicit  solutions 
as  comparison  functions.  This  will  be  done  briefly  in  Section  3  using  ideas 
of  V&zquez  [VI, 2].  The  principal  result  is 

Theorem  1.2.  Let  v  be  the  solution  of  (P),  where  the  smooth  datum  g 
satisfies  assumptions  (g),  xg(x)  >  0  for  x  near  zero,  and  there  exists 
Y  >  0  such  that 

-lim  g(x) |x|  -  a  >  0 

(1.7) 


x+0 


lim  G(x)x”^+1^  *  b  >  0  , 
+ 


x+0 


where  G(x)  *  f  g(£)d£,  x  >  0.  Let  p  :=  b(Y  +  1)/a.  Then 
0 

s(t)  -  -KtV2(1  +  o(1))  (t  ♦  0+)  , 
where  <  *  *(p)  >  0  is  the  uniquely  defined  inverse  of  the  function 


(1.8) 


VM 


As  already  remarked  in  the  special  case  y  =  1  the  integral  equation  approach 
yields  the  stronger  result  (1.1).  The  proof  of  Theorem  1.2  is  sketched  in 
Section  3. 

Proposition  1.1  and  the  comparison  method  used  in  the  proof  of  Theorem 
1.2  can  be  used  to  discuss  the  regularity  of  the  free  boundary  near  t  =  0 
for  initial  data  more  general  than  those  satisfying  (1.7).  For  example,  if. 


in  place  of  (D), 


ax  ,  x  >  0 

-|x|8,  x  <  0  , 


(1.9) 


where  a,  a  >  0,  8  >  0,  and  a  >  8,  then  there  exist  constants  c1rc2  >  0 
which  depend  on  a, 8  such  that 


c^a®+1t^  <  -s(t)  <  c2a^+1t^  , 


(1.10) 


Q  ^  1 

where  A  *  o7~ft  *  for  t  >  0  sufficiently  small.  On  the  other  hand,  if 


8  >  a  in  (1.9) , 


iim  .;2  ~s (-V. -  -  (2(8  -  a)) 1/2 

tV  t 1/2 1  log  t|1/2 


If  8  ■  a  in  (1.9),  A  ■  and  the  situation  is  covered  by  Proposition  1.1. 

The  constitutive  function  $(v)  *  max(v,0)  in  (P)  is  piecewise  linear. 
However,  self  similar  solutions  of  (P,D)  exist  for  more  general  constitutive 
functions  t,  but  they  cannot  in  general  be  found  explicitly.  We  expect  to 
investigate  the  regularity  of  the  free  boundary  near  t  *■  0  for  such  problems 
in  future  work.  Ibis  question  is  not  discussed  in  [BMPe]  where  other 
qualitative  aspects  of  several  such  problems  are  investigated. 


■W-T‘ 
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2.  SELF  SIMILAR  SOLUTIONS.  In  this  section  we  consider  the  Cauchy 
problem  (P)  for  the  model  datum  (D).  Let  v  =  v(x,t ;Y»P+»P_)  be  the  solution 
of  (P),  (D),  and  define  the  transformation  v  +  Tv  by 

Tv( x, t )  =  pv(Xx,X2t)  ,  (2.1) 

where  p,  X  >  0  are  constants.  T  transforms  solutions  into  solutions  and  in 
particular,  it  is  clear  that  (p_)v(x,t;y,p, 1 ) ,  where  p  *  p+/p_,  is  a 
solution  of  (P),  (D).  Denoting  v(x,t,y>p,1)  :m  v(x,tjy,p),  uniqueness 
(Theorem  3.2)  implies  that 

v(x,t»Y/P+,p_)  *  (p_)v(x,tfY,p)  .  (2.2) 

Since  Tv(x,0jy,p)  =  v(x,0;Y,p)  if  and  only  if  pX^  =  1,  uniqueness  of 
solutions  of  (P),  (D)  further  implies  that 

v(x,t;Y,p)  =  (X-Y)v(Xx,X2t>Y,p)  (2.3) 

for  any  X  >  0.  Fixing  the  point  (x,t)  and  choosing  X  =  t-1^  in  (2.3) 
implies  that  v  is  the  self  similar  solution 

v(x,t»Y,p)  s*  tT/2\|»(xt"1/2)  (2.4) 

of  the  Cauchy  problem  (P),  (D)  for  any  y  >  0.  By  the  comparison  method 
(Section  3)  the  free  boundary  s  of  (P),  (D)  is  monotone  decreasing  and 
s' (t)  <  0  for  t  >  0,  s(0)  =  0.  Moreover,  by  the  equivalence  of  problems 
(P),  (D)  and  (f),  (D)  the  free  boundary  s  ■  {(x,t)  :  xt"1//2  *  -k,  t  >  0} , 
where  <  >  0  is  the  constant  uniquely  determined  in  Proposition  1.1,  and  by 
Proposition  3.3 

v(x,t;Y»p)  >0,  x  >  s(t),  t  >  0  . 


Proof  of  Proposition  1.1.  Substituting  v  defined  by  (2.4)  in  (P)  one 
sees  that  lji  must  satisfy  the  linear  differential  equation  (1.3),  and 
conditions  (1.4)  and  (1.5). 

Equation  (1.3)  can  be  solved  explicitly.  Put  x  -  £//2  and 
w(£)  •  <|>(x).  Then  (1.3)  becomes 

w"(x)  +  xw'(x)  -  ywlx)  ■  0  .  (2.5) 

Setting  w(x)  «:  y(x)exp(-x2/4)  we  obtain 

2 

y"(x)  -  (y  +  Y  +  f-)y(x)  -  0  .  (2.6) 

This  differential  equation  has  the  general  solution  [B,  p.  116-117] 

y(x)  -  +  b2DY*iX*  <  *  <  **  7  >  0)  • 

where  Dy ( « )  is  the  parabolic  cylinder  function  of  index  v.  Thus  the 
general  solution  of  (1.3)  is 

♦t'C)  -  +  b,,DT(^)]exp(^-)  (2.7) 

for  —  <  £  <  •  and  y  >  0.  To  impose  the  initial  conditions  (1.4)  we  need 
the  formulae  [B,  p.  119] 


Then  the  initial  conditions  (1.4)  yield  the  pair  of  equations 


-7 


JH» 


•jBij 


J. 


Because  (2.6)  is  of  self-adjoint  form  the  Wronskian  of  D^(»),  D_^  ^(»)  is 


constant , 


_1(*),D  (•))  2  -iexp[ (-  *  1  )iri ]  . 


(p_,  J+1Dv(-  ^)exp(ll) 


b  1  <  tc ) 


i/2  exp[(l-t-l)wi] 


(2.9) 


b2(x) 


i/2  exp[  (Y  *  ^Tti] 


and  (2.7)  with  bi,b2  given  by  (2.9)  is  the  solution  of  (1.3)  satisfying  the 


initial  conditions  (1.4).  To  compute  the  limit  in  (1.5)  we  use  [B,  p.  122): 


D  (*)  -  zVexp(-  y-)[1  +  0(|z|“2)]  as  |z|  +  »  , 


(2.10) 


which  is  valid  for  -  <  arg  z  <  —■  i  Thus  for  £  e  R,  Y  >  0, 


=  exp(-  f") fe  Y  M  +  Oder2)],  C 

Y  W2  8  /2 


(2.11) 


■>  -  «xp(^)fe''  exp(iH)  „  .  oiler2)),  5 

T  So  So  * 


+«  . 


Substitution  of  (2.11)  and  (2.9)  into  the  general  solution  (2.7)  yields 


♦  (C)  “  b  (K)exp(i^-)(1-)  [1  +  0(1)],  ?■*■+•. 

i  /2 


(2.12) 


i*  J-  .*  , 


*.*  *.•  V  v*  •  V  *>  v  V  *>  V  V  W  V  ’  **  V  *>  *•'  •  ”  -  v  v  •  w  •  -  *  k>w>  / 


From  formula  (2.9)  we  see  that 


b2<k)exp(^)  =  ^  J  +  1D  (-  ^)exp(|— )  .  (2.13) 

/2  '  /2 

Imposing  the  asymptotic  condition  (1.5)  and  using  (2.12),  (2.13)  we  finally 
obtain 

lim  1LLL  =  p  =  p  £-)  D  (-  -^)exp(-~)  (2.14) 

?♦+-  5Y  Jl  ^  1  /2 

which  is  equation  (1.6). 

To  complete  the  proof  of  the  Proposition  we  have  to  show  that  given  any 
p  >  0,  (E)  is  uniquely  solvable  for  ic.  This  is  a  consequence  of  uniqueness 

(Theorem  3.2).  However,  it  can  also  be  seen  directly  as  follows.  From  [BO, 
p.  573] 


D  A-—) 


ft 


2n 


7+1 

z  2  r(i  +  I) 


n=0 


(2n)  J 


py 
/  2 


2n 


+ 


/ir 


00 


I 

n=0 


a2n+1 
(2n  +  1) ! 


2n+1 

9 


(2.15) 


is  an  analytic  function  of  <,  -»<*:<<»,  y>-1;ao  =  a.)  =  1, 

an+2  =  (7  +  2-)an  +  f  (n  -  1)an_2,  and  D_y_,(0 

K 

a 

increasing  function  of  <  for  0  <  x  <  »,  and  by  (1.6)  so  is  p(tc). 

Moreover  p(0)  =  0.  This  completes  the  proof  of  Proposition  1.1. 

We  conclude  this  section  by  establishing  two  useful  asymptotic  estimates 


the  coefficients  a  are  positive,  D  .  f 

v  -Y-1'- 


2 

is  a 


(7+1 )/2 


r(i  + 

positive. 


B 


Since 


strictly 


From  (2.15)  and  ( E ) 


Moreover,  [BO,  p.  574] 


%-i(-  *  r TTTyT  ■ 


and  therefore,  from  (E) 


rz —  2y+1  z 

P(IC>  "  rrrrTT  (y  ®*p(t-)  (<■*•+»)  .  (2.17) 

'  V2 


3.  EXISTENCE  UNIQUENESS  AND  COMPARISON  OF  SOLUTIONS  OF  (P).  In  this 

section  we  state  general  existence  and  uniqueness  results  for  (P),  briefly 

discuss  the  comparison  method,  and  sketch  the  proof  of  Theorem  1.2. 

Consider  the  Cauchy  problem  (P)  in  f^-Mx  [0,T)  where  T  >  0,  and 

00 

assume  that  the  datum  g  e  and  satisfies 

g(x)  <  c  exp(hx2),  x  e  R  ,  (3.1) 

for  some  constants  c,  h  >  0.  Define  g“(x)  :=  max( -g(x) , 0 ) .  For  existence 
of  solutions  of  (P)  one  has  the  following  result  proved  in  V&zquez  [VI]. 
Theorem  3.1.  There  exists  T  ■  T(g)  >  (4irh)  1  and  a  function 
v  e  C(  [0,T]  :  L^q^R)  O  L1oc(Qt))  with  the  following  properties: 

(i)  v(»,t)  +  g  _in  l|oc(R)  as  t  -*■  0+. 

(ii)  vt  -  ♦(v)^  =  0  iji  V  (&!>)• 


(iii)  For  every  t <  <  T  there  exists  constants  h«  >  h  and  c,  >  0  such 


Uniqueness  of  solutions  of  (P)  can  be  established  in  a  larger  class  of 
functions  which  includes  the  class  of  solutions  in  Theorem  3.1;  the  proof  is 
given  in  VSzquez  [VI]. 

Theorem  3.2.  Let  u,v  c  C([0,T]j  n  L1oc(Qt))  satisfy  the  following 

properties : 

(i)  ^  -  ♦<u)xx  *  vt  -  $(v)xx  in  t?'(QT); 

(ii)  u(*,t)  -  v(»,t)  *  0  in  I»ioc(*)  t  ♦  0+; 

(iii)  for  every  t-j  <  T  there  exist  constants  c<|  >  0,  hj  >  0  such  that 


u(x,t) ,v(x,t)  <  Cjexp(h^x2)  a.e.  _in  Rx  (0,t1) 


(3.3) 


(iv)  u"(x,t),v”(x,t)  €  l1oc<Qt)‘ 


Then  u(x,t)  *  v(x,t)  a.e.  ^n^  Qp. 

As  a  consequence  of  Theorems  3.1  and  3.2  the  Cauchy  problem  (P)  has  a 
unique  solution  in  Qbp  for  sosie  T  >  0  if  the  datum  g  satisfies  (3.1).  If 
the  datum  g  has  exponential  growth  as  |x|  ♦  •,  the  solution  v  of  (P) 
will  become  infinite  in  finite  time.  This  is  evident  because  nonnegative 
solutions  of  the  heat  equation  vt  -  vxx  have  this  property  for  exponentially 
growing  data.  On  the  other  hand,  if  the  datum  g  has  at  most  polynomial 
growth  at  infinity,  then  by  Theorem  3. 1  the  solution  v  of  (P)  exists  on 
0  <  t  <  “». 

Next,  we  turn  to  the  comparison  method  in  the  form  it  is  needed  in  the 
proof  of  Theorem  1.2,  see  VSzquez  [VI].  To  state  it  we  define 

x 

G(x)  -  /  g(5M5»  x  e  R  , 


V(x,t)  -  /  v(£,t)d£,  x  €  R,  t  >  0 


Wm 


Proposition  3.3.  Let  Vi.Vi  be  two  solutions  of  (P)  corresponding  to  data 


g.|,g2  respectively,  each  satisfying  assumptions  (g),  and  9ft 92  e  L 
G.j(x)  >  G2(x)  for  every  x  e  R,  then  V.j(x,t)  >  V2(x,t)  for  every  x  e  R, 
t  >  0. 

Remark.  Proposition  3.3  holds  for  solutions  v  of  (P)  whenever  t  is  a 

continuous  nondecreasing  function:  R  R,  and  $(0)  =  0.  Observe  that  V 

satisfies  the  P.D.E.  V*.  =  A(V„)  and  Proposition  3.3  is  the  maximum 

*  x 

principle  for  this  equation. 

Sketch  of  Proof  of  Theorem  1.2.  We  assume  that  assumptions  (g)  are 

satisfied.  Since  the  interest  is  in  the  regularity  of  the  free  boundary 

near  t  ■  0,  we  assume  without  loss  of  generality  that  g  has  at  most 

polynomial  growth  at  infinity.  Then  the  Cauchy  problem  (P),  (g)  exhibits  a 

unique  free  boundary  x  *  s(t)  defined  on  0  <  t  <  »,  s( 0)  =  0,  s  e  C^(0,«»), 

and  s'(t)  <  0  for  t  >  0.  In  what  follows  we  denote  by  G(x)  the  primitive 
x 

/  g(5)d£.  An  essential  ingredient  in  the  proof  is  the  following  refinement 

0 

of  Proposition  3.3.  An  analogous  result  for  the  porous  medium  equation  has 
been  proved  by  V&zquez  [V2] »  the  proof  is  similar. 

Lemma  3.4.  Let  v,.v,  be  two  solutions  of  (P)  corresponding  to  initial 


data  g<|,g2  respectively,  satisfying  assumptions  (g).  Let  s^,s2  be  the 
corresponding  free  boundaries,  s^(0)  *  s2(0)  =  0.  Assume  that 


for  small  x  <  0,  and  that 


0  >  gj(x)  >  g2(x) 


G«(x)  >  Go(x) 


for  every  x  on  some  interval  (0,a),  o  >  0,  G^ta)  >  0  and  Gj  £  G2 


(0,a).  Then 


s^(t)  <  s2(t) 


for  t  >  0  sufficiently  small 


V 


To  complete  the  proof  of  Theorem  1.2  we  apply  Lemma  3.4.  Let 

the  solution  of  (P),  (g)  where  the  initial  datum  g  satisfies  assumptions  (g) 

and  (1.7).  Let  V2  be  the  self  similar  solution  of  Proposition  1.1 

corresponding  to  the  model  datum  (D)  with  p+=  b(  y  +  1)  +  e ,  p_=a  -  t , 

where  0  <  e  <  a.  Then  G2<x)  >  Gj(x)  for  small  x  >  0,  t  G 2»  and 

0  >  g2(x)  >  g<|(x)  for  x  <  0.  By  Lemma  3.4  s2(t)  <  s^Ct)  for  small  t  >  0. 

Thus  identifying  Sj  with  s  and  using  the  result  of  Proposition  1.1, 

liminf  (-s(t)t1/2)  >  <(b(Y  +  1}-+  C) 
t.0+ 

for  every  0  <  e  <  a.  Since  ic  is  a  continuous  function  of  p+/p_ ,  we  let 
e  ♦  0+  to  obtain 

liminf  -t-?1)  .  (3.7) 

t.0+ 

By  a  similar  argument,  letting  be  the  self  similar  solution  of  (p),  (D) 

with  p+  and  p_  defined  as  above,  and  letting  v2  be  the  solution  of  (P), 
(D)  we  obtain 

limsup  (-s(t)t1^2)  <  .  (3.8) 

t-0+ 

The  result  (1.8)  is  a  consequence  of  (3.7),  (3.8).  This  completes  the  sketch 
of  the  proof . 
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